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A disorder-dependent Gaussian variational approach is applied to the 
d-dimensional ferromagnetic XY model in a random field. The randomness 
yields a non extensive contribution to the variational free energy, implying a 
random mass term in correlation functions. The Imry-Ma low temperature 
result, concerning the existence {d > 4) or absence [d < 4) of long-range order 
is obtained in a transparent way. The physical picture which emerges below 
d = 4 is that of a marginally stable mixture of domains. We also calculate 
within this variational scheme, disorder dependent correlation functions, as 
well as the probability distribution of the Imry-Ma domain size. 
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I. INTRODUCTION 



The effect of quenclied disorder on systems with continuous symmetry has recently at- 
tracted a lot of theoretical as well as experimental interest. Experimental realizations in- 
clude, among others, arrays of fiux-lines in type II disordered superconductors, crystalline 
surfaces with a disordered substrate, spin or charge- density waves subject to random pinning, 
etc... Most theoretical studies of this model have focused on the random-field vortex- free 
XY case, which is equivalent to the random phase sine-Gordon model. P|-P^ 

Recently, a disorder dependent variational approach has been proposed for this prob- 
lem |]18|. In this approach, the disorder enters only through a unique variable u = 



J dxcos(2TTd{x)), where d{x) is a random phase. This random variable has a Gaussian 
distribution, and for u < 0, one recovers the results of the replica Gaussian variational 
principle, with Parisi symmetry breaking scheme 0. 

One major advantage of this approach is that it is genuinely variational, thus providing 
a true upper bound to the free energy of the system, unlike replica based methods, which 
are plagued by the n = limit. 

In this paper, we shall use this variational approach in the framework of the (full) XY 



model in a random magnetic field, for which few results are available [T^JSy] . Its Hamiltonian 
reads: 

H = — J ^ cos{9f' — Ofti) — ^ h^cosiOfT — 0^) (1) 

<r,r'> r 

where denotes the phase angle of the XY spin S^^ J is the ferromagnetic nearest-neighbor 
coupling constant, and {h^^ (pi) are the polar coordinates of the random field hf^. The prob- 
abihty distribution of the (site uncorrelated) random field is Gaussian, given by: 

where N = L'^ is the total number of sites on a dimensional hypercubic lattice of linear 
size L and lattice spacing a; h is the variance of the random field. 



The physics of the Hamiltonian (P is beheved to be well captured by the Imry-Ma 



argument, which we briefly summarize |T^. Consider, at low temperature, a magnetized 
domain of size L, with magnetization in the direction. To study the stability of such a 
region, we imagine that we construct inside a subdomain of size ^, in which the magnetization 
is aligned with the average local magnetic fleld. The energetic balance reads: 

~ - (3) 

where the flrst term represents the spin-wave distortion energy, whereas the second represents 



the magnetic energy gain due to the random fleld. Following we conclude that the lower 
critical dimension for the system is dc = 4. Above dc, the energetic cost is prohibitively high, 
so that such domains of size ^ cannot exist: ferromagnetic long range order is stable. On 
the contrary, for d < dc-, the whole system will break into sub domains of smaller size, so 
as to maximize its magnetic energy gain. Consequently, no long range ferromagnetic order 
may exist below d^.. 

Usually, in quenched disordered systems, extensive thermodynamical quantities, such as 
the free energy, are identifled to their average over the disorder This can be understood 
in the following way: one divides the macroscopic system into mesoscopic subsystems , each 
subsystem corresponding to a particular disorder conflguration. For short-range forces, the 
free energy is additive, and thus the total free energy is the sum of the free energies of the 
subsystems. This procedure clearly neglects all correlations or domain-wall energies between 
neighboring subsystems. In the present case, these non-extensive contributions are precisely 
of the same order of magnitude than the terms of equation (^. Therefore, in the following, 
we will not perform quenched averages, but rather keep the disorder variables throughout 
the calculations. 

The layout of this paper is the following. In section II, we deflne the variational Hamilto- 
nian, and calculate the corresponding variational free energy. The variational method yields 
two solutions. In this approximation, we flnd that the transition temperature is the same 
for both solutions, equal to that of the pure system. In section III, we discuss the issue of 



long-range order as a function of space dimension. We find long-range ferromagnetic order at 
low temperature in dimensions d > A. For lower dimensions, the discussion is postponed to 
the next section. In section IV, we discuss the non-extensive corrections to the free energy, 
and study the stability of the variational solutions. For dimensions 2 < < 4, we argue that 
the physical solution is the marginally stable one, shedding light on the Imry-Ma domain 
picture. In the conclusion, we calculate the probability distribution for the Imry-Ma length, 
and discuss the issue of correlation functions. 



II. THE VARIATIONAL FREE ENERGY 

We consider Hamiltonian (|l|) and its associated variational Gaussian companion: 

/^^O = I J2i(^r - Oo) 9if- f'r'ie,, - 60) (4) 

where we have restricted the variational kernel g to be translationally invariant, the direction 
of magnetization 60 is space independent and j3 = 1/T is the inverse temperature. 
The true free energy F satisfies the bound: 

F < $({/!,-}) = Fo+<n~no >o (5) 

Using equations we obtain the disorder dependent variational free energy: 

1 - 7 

Hihr}) = -^T. lii^(^) - ijL'e-^^'^ J2 e^^^-^ - e'^^/^ ^ /i.^cos(0.^ - ^0) (6) 
^ k "=i r 

where g{k) is the Fourier transform of g{r), {cq,}, a = 1, ...,2d denote the lattice unit vectors 

and {hf, 0^} are the polar coordinates of the random field hf at site f. 

The variational equation with respect to 6*0 reads: 

e-5(o)/2 j2 hr-sm{(t),. - 0o) = (7) 

r 

If we define as the polar angle of the total magnetic field H = J2rhf, we see that the 
solutions to (|^ are 6*0 = 4>, corresponding to a magnetization aligned with the total magnetic 
field H, and 6q = (p + tc corresponding to a magnetization opposite to H. 

4 



The minimization of ^{{h^}) with respect to g{k) then yields, for the two solutions: 



1 2d 

— ^ = J5](1 -e''=")e-(^'(°)-^("")) ±e-55(o)L-'^^/i,-cos(0,--0) (8) 
pg{k) a=i 

the + (resp. — ) sign corresponding to 6*0 = (resp. 6*0 = + vr). For the sake of simplicity, 
we shall stick below to a unique notation g(k) for the two solutions. In the following, we 
shall refer to these solutions as the (+) solution (for 9o = (p) and the (— ) solution (for 
^0 = + vr). 

Defining the components = X^r cos(0,^) , hy = J2fhfsin{4>r) , we see that: 

J2hr-cos{^r^- (j)) = H (9) 

r 

— > 

The total magnetic field if is a random variable, and the distribution of its modulus is given 
by the central limit theorem. Defining the positive random variable u by: 

H = L^'"^ hu (10) 

we have: 

P(m) = Mexp(-MV2) (11) 

Each disorder configuration is thus specified by a single positive random variable u. The 
u-dependent free energy per site (the upper sign corresponding to the (+) solution, and the 
lower sign to the (— ) solution). 

m = ^ = / - E e^(-) T e-^^^y-L-^'-hu (12) 

where we have replaced by L'^ j (see below). 

Using the symmetry of the problem, we set 71 = g{Qi) — g{ea) independent of a. Summing 
over a, we get: 

2d 2d Jdu 

dli = T.i9{0) - g{ea)) = Y: TTTT-Jim - e^'-) (13) 

Inserting equation (^) in equation ([TSp, we obtain the (disorder dependent) critical temper- 
ature through: 



2rfJe-^i7i = Tt g{0)e-y^^^L-'^/^hu (14) 

For large system size (L — oo), we see that the second term of the r.h.s is always negligeable, 
and is therefore identical for the two solutions. The critical temperature is the same as in 
the pure system T^^^"^ = where e = 2.718... The quantity 71 varies from (at T = 0) 
to 1 (at TP"''^). 

As is clear from equation (^, the form of g{k) and therefore of 0(m) is very different for 
the two solutions. For the (+) case, the replacement of ^'^ I (Cj^ straightforward, 

whereas the (— ) case requires a more careful treatment, due to the existence of poles in g{k). 



III. LONG RANGE ORDER AND THE IMRY-MA ARGUMENT 

A. The (+) solution 

In this case, equation (|^) shows that: 

{2txY I ^2d^^(i - eifc.)e-(9(o)-9(e.)) + e-'^^io) ^^^^ 

implying (for d> 2) the existence of long range order both in the angular and spin variables. 
Indeed, we have: 

<(^^r-^^o)'>=^7(0) (16) 

< >= COS0 exp(-^) (17) 

and (?(0) is finite (infrared convergent) for d > 2. 

For d = 2, one finds, as in the pure case, algebraic order, but with different exponents. 
For example, the total magnetization reads: 

= < >= e-^ ~ (^^^ COS0 (18) 

with ^ = ^^e""^^, the pure system behaving as 



M^ure - [-l) (19) 

Note the appearance of the {u dependent) Imry Ma length scale ^ in equation |l^. We do 
not investigate further the case (+), since it will be shown, in section |^ that the physically 
relevant case (2 < (i < 4) is the (— ) one. 



B. The (-) solution 

In this case, one has to be more careful in taking the continuous limit in (p!2D, since the 
poles may yield a finite contribution to thermodynamic quantities (in close analogy to the 
Bose condensation). For instance, we may write (for small wave vectors): 

/3JL'^^;(0)=E-^ (20) 

where 

^ = ^ (21) 

In equation (|20|), there may arise a singular contribution from the 2d smallest wavevectors 
k = ^(1, 0, 0, ...) (plus permutations). One therefore has 

PJL''g{Q) = +L' T^,^ (22) 



A detailed discussion of equation (p^ requires separate treatments for (i > 4, for 2 < c/ < 4 
and for d = 2, 4. 



1. Dimension d > 4 

According to (pT]), we have fi^ << L~^. We may thus neglect the /i^ contribution in the 
denominators of (P^, yielding a finite magnetization: 

= — X] < 55 >= -e'Ti^) COS0 (23) 



where A 



ell 



Ftt^ T{d/2)a<l 



2. Dimension 2 < d < 4 



Since the integral in (^) is infrared (IR) convergent, it may be neglected compared to 
the pole contribution. In this case, one finds: 

(iJL''g{^) ^ (24) 



This last equation shows that (when L oo), /x is very close to 2n/L. Using equation (pi|), 
we get: 

(7(0) = (4-rf)ln-^ (25) 



where 



Since g{0) diverges for large L, we conclude that there is no long range order in this case. 

3. Dimensions d = 2 and d = 4 

For d = 2, one may easily see that both terms in the r.h.s. of equation (^) are of the 
same order of magnitude. Indeed, the integral has a logarithmic divergence, and performing 
the calculation, we obtain the same result as above, namely: 

^?(0) = 21nf-^^ (27) 



with 



For c? = 4, one gets: 



implying disorder induced fluctuations in the magnetization. 



«">=H^ '''' 



m-M^-^\ (29) 



IV. MARGINAL STABILITY AND FINITE SIZE CORRECTIONS 



Whenever there are several solutions to the variational equations, one ought to pick the 
one with the lowest free energy. In our case, it is easily seen that the extensive part is the 
same for both solutions, and that the difference shows up only in non-extensive corrections. 

A. Dimension 2 < d < 4 

To leading order, the variational free energy reads, 

= FqL'^ - A+ hu L'^/^ (30) 



= FqL'^ + (4 - rf) A_ Je~^i L'^-'^ In ( — ^ ) 

\idiu)) 



(31) 



Fo = -i- In f - dJe-^' + / ln(fc2) (32) 



where 



2(5 \(3J) 2(3 J {27tY 

and and are positive constants, independent of L, depending on geometrical factors 
as well as on the temperature of the system. 

The above selection criterion holds only if the two minima are not degenerate. In the 
opposite case, the selection rule is provided for by the fluctuations since they may contribute 
to the extensive part of the free energy. In other words, one has to check for the stability of 
the solutions, with respect to variations of g{k) and local variations 6ff of ^^o- 

We consider first the 5^ fluctuation contribution to the free energy and get: 



/ J 2d -, 

= E 46""^ E i^r - WJ' + 2 cos(# - ^o) 5' 

r \ a=l 



(33) 



The positivity of AF is determined by the spectrum of the kernel in (pSj). This kernel is 
analogous to that of a tight-binding model in a random potential: 

-Za,-,-, + V^^ = (34) 



where 

V,. = e-f (°)/2e7i cos{(f)r- - Oq) (35) 
It is easy to see that is a random Gaussian distributed potential with: 

V^=e-9(^y^^hL~''/^ -^0 (36) 

V^ = e-^'(°)/i25(f-0 (37) 

where the bar stands for an average over the random field. 

Using standard perturbation theory to second order, the ground state energy of (p^ can 
be expanded around the zero energy mode k = as: 

=< 0|y|0 > - V \<f'\^\^>\' (38) 

where \k > denote the normalized plane waves. 

Since is a random potential of mean given by (^61), the first term of the r.h.s. of (|38D 
is of order V ~ L~°'/^. The second order term is negative and finite, and thus the energy 
Eq is negative, leading to a fiuctuation induced instability of the finite magnetization (+) 
solution. 

This instability is in agreement with the Imry-Ma picture, which shows that for d < 4, 
a uniformly magnetized system is unstable with respect to domain formation. 

On the contrary, for the zero-magnetization (— ) solution, the disordered potential van- 
ishes. The spectrum is that of a free particle, implying the existence of zero energy modes. 

Furthermore, one may check that the (— ) solution is stable with respect to variations of 
g{k), since its stability properties with respect to these variations are the same as the pure 
system. 

Putting all these results together, we may conclude that the (— ) solution is marginally 
stable, and thus appears as the physical solution. 
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B. Dimension d > 4 



It is easily seen that both solutions having finite magnetization, (aligned or opposite to 
the field), they are unstable. This is again in accord with the Imry-Ma argument, which 
shows that for d > 4, the magnetization (i.e. ^o) is not determined by the random field, but 
rather by a standard infinitesimal uniform field. This in turn ensures the stability of the 
uniform magnetisation with respect to small fiuctuations. 



C. Dimensions d = 2 and d = 4 



According to section |^, for d = 2, both solutions are marginally stable. This is a 
borderline case, due to the existence of a Kosterlitz-Thouless transition in the pure system. 



For = 4, it can be seen that the random potential of equation (34) is not Gaussian. 
This dimension is likely to be a special dimension for the localization problem. This fact is 
also present in the framework of the d ^ d — 2 dimensional reduction theory [E2 , 23 . 



V. CONCLUSION 

We have presented a disorder dependent variational method for the full XY model in a 
random field. We recover the Imry-Ma results concerning the existence {d > 4) or absence 
(2 < (i < 4) of long range order. In particular, we find in the latter case, a variational solu- 
tion which is marginally stable with respect to local magnetization rearrangements. Using 



eqs.(|rT|) and (|26|) , this solution can be further characterized by the probability distribution 



of the Imry-Ma domain length ^d, which reads 



, (4-d) 5 , B , 



with B = ( 47r je Ti Eq.([39|) shows in particular, that the Imry-Ma domains are charac- 
terized by multiple length scales; for d=3, the average domain size is found to be divergent, 
whereas the most probable domain size is of order B. 
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In a similar way, one may deduce the spin-spin correlation function. We have 

K{r) =< SoSr- >= e-(^(°)-^("^)) (40) 
Considering only the case 2 < c? < 4, we have 



)ln- 



K{r}=e V ' ' ^'^'"V (41) 

where 

^^) = X:(l-cos^) (42) 

and Xa denote the coordinates of r. The disorder dependence of the correlation function 
is contained in the Imry-Ma length scale ^^(m): one may therefore get its probability dis- 
tribution. Here we just point out two limiting cases for the disorder dependent spin-spin 
correlation function. At small distances {xa << L), K{f) is Gaussian. At large distances 
{xa ~ L), K{r) decreases like a power law 

K{r) ^ (^)-^ (43) 

where w is the value of w{^) for Xa ~ L. This behaviour is in broad agreement with the 
results of (i) a real space renormalization group (ii) a replica variational calculation 0, 
both calculations pertaining to a vortex free model. The full XY model in a random field has 
been recently studied by Monte Carlo calculations for d = 2,3, indicating (for d = 3) 
the possibility of a phase transition to a pinned vortex free phase. Our variational approach 
emphasizes the existence of a probability distribution for the Imry-Ma domain size, which 
renders the comparison with this work rather delicate. Finally, our results should be of 



interest in other Imry-Ma like situations p^j25|] 
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